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UNIFYING THE BGM AND SABR MODELS: 
A SHORT RIDE IN HYPERBOLIC GEOMETRY 

PIERRE HENRY-LABORDERE 



Abstract. In this short note, using our geometric method introduced in a previous paper 1121 
and initiated by [1], we derive an asymptotic swaption implied volatility at the first-order for a 
general stochastic volatility Libor Market Model. This formula is useful to quickly calibrate a 
model to a full swaption matrix. We apply this formula to a specific model where the forward 
rates are assumed to follow a multi-dimensional CEV process correlated to a SABR process. 
For a caplet, this model degenerates to the classical SABR model and our asymptotic swaption 
implied volatility reduces naturally to the Hagan-al formula fill . The geometry underlying this 
model is the hyperbolic manifold H" +1 with n the number of Libor forward rates. 



1. Introduction 

The BGM model jHlE] has recently been the focus of much attention as it gives a theoretical 
justification for pricing caps-floors using the classical Black-Scholes formula. The basic (physical) 
random variables arc given by the Libor forward rates which are assumed to follow a correlated 
log-normal process. As the forward swap rate model implied by the BGM is quite complicated 
(the swap forward rate is not log- normally distributed), the calibration to the swaption matrix 
is difficult. Asymptotic swaption implied volatility (at the zero-order) were initially derived by 
Rcbonato and Hull- White |H| for the (log-normal) BGM model. Such formula has been obtained 
by assuming that the ratio of a forward Libor rate over the swap rate and the derivative of the swap 
rate according to a forward Libor rate are almost constant (and therefore equal to their values at 
the spot). 

Despite its great success, the BGM model presents the same drawbacks as the classical Black- 
Scholes theory: As the forward rates follow a correlated log-normal process, the model is not 
able to calibrate the full swaption matrix in/out-the money (in particular the caplcts) and give 
a good dynamics to the Libor rates. The incorporation of a swaption smile can be obtained by 
introducing more elaborated models which should be flexible enough to calibrate caplets and a grid 
of swaption volatilities (not necessary at the money) across all swaption expiries and underlying 
swap maturities. One property that these models must still share is their ability to quickly calibrate 
the swaption matrix without using complicated numerical routines such as Monte-Carlo simulation 
which are usually noisy and time-consuming. In this context, Andersen- Andreasen introduced 
the CEV Libor Market Model (LMM) P which assumes that each forward rate follows a CEV 
process, and showed how to obtain asymptotic swaption smile. Their method is still based on the 
Rcbonato "freezing" argument which is not completely mathematically justified. Recently, for this 
specific model, Kawai found a better asymptotic formula using the Wiener chaos expansion |16| . 
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Although giving more flexibility than the BGM model, the CEV LMM model is still not able to 
calibrate the swaption matrix for in/out strike and in this context, we are naturally led to use 
stochastic volatility LMM. The literature on this subject is not particularly large. Andersen-al 
introduced a LMM where the Libors follow a multi-dimensional correlated CEV process coupled 
(but uncorrelated) to a Heston model [2| [3| and recently V. Piterbarg modifies this model to 
incorporate term structure |f 7| . Using an averaging principle, V. Piterbarg derives an asymptotic 
volatility. Note that as these models are uncorrelated to the stochastic volatility, the swaption fair 
value is simply given by the fair price in the case of a local volatility model conditional to the 
stochastic volatility process as explained by the Hull- White decomposition |H3]- An asymptotic 
expression can then be generated by approximating the moments of the volatility process |2]. 

For pricing exotic options (such as bermudan swaptions for example) , it is simpler or more natural 
to model directly the forward swap rate with a stochastic volatility process. For example, the 
SABR model |T1] was introduced to fulfill his goal. An asymptotic swaption smile formula (at the 
first-order) was derived for this specific model and help to calibrate quickly the model to liquid 
market data. In this context, it is natural to try to reconcilate/unify the two benchmark models, 
the BGM and SABR models. We therefore introduce a LMM where the forward rates follow a 
multi-dimensional CEV process (with one beta for each forward) correlated to a SABR model. 
As it is the case for the SABR model, we impose that the libors are correlated to the unique 
volatility and it is therefore not possible to follow the Andersen-al [3 method (i.e. the Hull- White 
decomposition) to derive an asymptotic swaption smile. 

In this paper, we pursue our previous work on the application of the heat kernel expansion on 
a Ricmannian manifold endowed with an Abclian connection |12{ to derive an asymptotic smile 
formula for a swaption. The plan of this paper is as follows: In the first part, we will recall some 
definitions and present a list of recent Libor Market Models. In the second part, we apply this heat 
kernel expansion to derive an asymptotic swaption smile formula at the first-order valid for any 
LLM. In the third part, we present our stochastic LMM and apply this general formula. We will 
prove that the geometry underlying this model is the hyperbolic manifold H™ +1 . Some important 
properties of this space are then presented. Furthermore, we show that the "freezing" argument 
is no longer valid when we try to price a swaption in/out the money: The libors should in fact be 
frozen to the saddle-point (constrained on a particular hypcrplane) which minimizes the geodesic 
distance on H n+1 . 



2. Libor Market Model 

A swaption gives the right, but not the obligation, to enter into an interest rate swap at a 
pre-determined rate on an agreed future date [7j- The maturity date for the swaption is noted T a 
and Tj3 is the expiry for the forward swap rate s a p given by 

(2.i) SaP (t)- 1 - n^ tt+1 ti+tW 



zCi= Q +i T Yij 



=a+l (1+TFjit)) 



with r the tenor and Fk(t) = F(t, Xfc_i, Tk) is the forward rate resetting at T^-i- 

As the product of the bond P(t,Tk) with the forward rates Fk(t) is a difference of two bonds 
with maturity XJt-i and Tk, -(P(t,Tk-i) — P(t,Tk), and therefore a traded asset, Fk is a (local) 
martingale under Q fe , the (forward) measure associated with the numeraire P(t,Tk). Therefore, 
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we assume the following driftlcss dynamics 

(2.2a) dF k (t) = a k {t)$ k {a,F k )dW k ,Vt<T k _ r ,k = !,-■■ ,n 

(2.2b) dW k dW t = Pkl (t)dt 

with the initial conditions a(t = 0) = a and F k (t = 0) = F®. 



Libor market model 


SDE 


BGM 


dF k = a k (t)F k dW k 


CEV 


dF k = a k (t)F^dW k 


Limited CEV 


dF k = a k {t)F k mm{F^,e^- 1 )dW k 
with e a small positive number 


Shifted log-normal 


dX k = a(t)X k dW k 
with F k = X k + a 


FL-SV 


dF k = a k (t)(PF k + (1 - p)Fj!)^dW k 
dv = X(v - X)dt + Uy/vdZ ; dW k dZ = 


FL-TSS 


dF k = a k (t){p k {t)F k + (1 - /3 k (t))F^dW k 
dv = X(v - X)dt + v^jdZ ; dW k dZ = 



Table 1. Examples of stochastic (or local) volatility Libor models. 



In order to achieve some flexibility, we assume that the (normal) local volatility & k (a, F k ) depends 
on a hidden Markov process a (to be specified later) representing a stochastic volatility. We 
therefore assume that all the forward rates are coupled with the same stochastic volatility a. 
(Table 1) presents a list of the different functional forms for <f> k used in the literature. The BGM, 
(limited) CEV and shifted log- normal models correspond to local volatility models (a = 1) and 
the others to stochastic volatility models with a unique stochastic volatility a driven by a Heston 
process. Note that the stochastic differential equation for the libors L k has been written in the 
forward measure Q k and the stochastic equation for a remains the same in the forward or forward 
swap rate measures as a is assumed to be uncorrelated with the Libor rates. This will not be the 
case in our LLM. 



3. Asymptotic Swaption Smile 



The forward swap rate satisfies the following driftless dynamics in the forward-swap measure Q) a/3 
(associated to the numeraire C a p(t) = Y^i= a +i T i-P(^^i)) 

13 q s 

(3.1) ds a/ 3 = ^2 -^r^k(t)4> k (a, F k )dZ k 

k=a+l 



A- 



with - S -^H I gfelk] _L £g=j TfcP ( t ' Tfc + 1 K 

Willi — (l+TjFj)^ P(t,Tp)-P(t,T a ) ~T C af3 I 



The local volatility associated to the forward swap rate (ds a p = c'[f c (sa0, i)dWi) is then by defini- 
tion 





(3.2a) «f) 2 (M) = E Q/3 [ £ pijWaiWvjWM^miM) 



ds a p ds a p , 
dFi dF 3 



\s a p = s\ 



(3.2b) Ptf( f M*)°-j(*) 

i,j—a-\-l 



ImPdaRidFi 
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with the submanifold B = {{Fi}i\s a/ 3 = s} and p = p(a,Fi,t\a,F®) the conditional probability 
satisfying the (backward) Kohnogorov equation associated to the SDE for the Libors and the 
volatility a in the forward swap measure Q"' 3 . An asymptotic expression in the short time limit 
for the local volatility a a 0(s,t) can be found in two steps: find an asymptotic expansion for the 
conditional probability p (in Q Q/3 ) and do the integration over B. 



3.1. First step: Heat Kernel expansion. As explained previously, the first step can be achieved 
using the heat kernel expansion. In that purpose, the Kolmogorov equation is rewritten as the heat 
kernel equation on a (n + l)-dimensional Ricmannian manifold A4 n+1 endowed with an Abclian 
connection as explained in |121 1131 HI] . Let's assume that our multi-dimensional stochastic equa- 
tions (in Q a/3 ) are written as 

(3.3) dx» = b»(x,t)dt + a"{x, t)dW" 

with dW^dW — p^v (t)dt (note that the indices 1, ■ • • , n (resp. n+ 1) correspond(s) to the forward 
F l (resp. a)). Then, the metric g^ v depends only on the diffusion terms and the connection 
on the drift terms as well 



p^(t) 



(3.4a) g^(x,t) = 2 ,\^' ( tx , H,t> = 1 ■ ' ' n + 1 , p^ = [p~^ 



J/"' 



n+l 



(3.4b) A*(x,t) = -(b^x^-^a^Mg^g^ix^))), (,= !■■ -n+l 



with g(x, t) — detlg^uix, t)]. In terms of these functions, the asymptotic solution to the Kolmogorov 
equation in the short-time limit is given by (x = (a, Fj), x = (a, F®)) 



(3.5) p(x, t\x°) = ^^JA(x,x°)V(x, x^e" 1 ^ 1 V a n (x, x°)t n , t -» 

(Ant) 2 z — ' 

• Here, a(x,x°) is the Synge world function equal to one half of the square of geodesic 
distance d(x,x°) between x and a; for the metric g^ v {x,t = 0). This distance is defined 
as the minimizer of 

f T dx^ dx v 

(3.6) d(x,x ) 2 = mine 9uv(x,t = 0) — — — dt 

J Q dt dt 

and t parameterizes the curve C{x, x°) joining x(t = 0) = x° and x(T) = x. 



A(x,x ) is the so-called Van Vleck-Morcttc determinant 



(3.7) A(x, x°) = g(x, 0)~i det(-^^l)g(x°, 0)"* 

with g(x, 0) = detlg^x, 0)] 
• V(x,x°) is the parallel transport of the Abelian connection along the geodesic C(x,x° 
from the point x° to x 



(3.8) V(x, x°) = e^*.* ) -*M(x.*=o)«to* 

• The ai(x,x°) (ao(x,x°) = 1) are smooth functions on M and depend on geometric invari- 
ants such as the scalar curvature R. More details can be found in |12j . 
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3.2. Second step: Saddle-point method. The integration over B is obtained by using a saddle- 
point method which consists in approximating at the first order the integral / / (x)e € ^^ dx in the 
limit e small by [8] 

f(x)e^dx ~ e «! f(x*)e^ {x *\l + ~{-^j-A a0 + (^fd^ + {d^^A^s 

e ZJ ZJ 8 

j2 A a pA lS A^)) 

with A a P = [dajjcf)]" 1 , dx = Y[?=i^ x i anc ^ x * the saddle-point (which minimizes <j>(x)). This 
expression can be obtained by developing <f>(x) and f(x) in series around x*. The quadratic part 
in <p(x) leads to a Gaussian integration over x which can be performed. 

d(x,x°) 2 

3.2.1. Saddle-point. As the conditional probability at the zero-order is proportional to e 4t , 
the saddle-point corresponds to the point x = (a, fi=i,... >n ) on the submanifold s a p = s which 
minimizes the geodesic distance d(x,x°) 0E] 

(3.10) (a*, {F*}) = (a, {Ft}) such as min a ^ Fi y\ Sal3=s [d(x, x ) 2 } 
Introducing a Lagrange multiplier, A, this is equivalent to 

(3.11) {a*, {F*}) = (a, {FJ) such as min a , {Fi} Jd 2 {x, x°) + X(s a0 (F) - s)] 

3.3. Asymptotic local volatility. Plugging our asymptotic expression for the conditional prob- 
ability H3.5fl into 13.9(1 . we finally obtain the local volatility at the first-order 

(af) 2 (s,t) = Y p ij (t)a i (t)a j (t)f ij (F*,a*)(l + 2t V A»» ( ( ^ 7 ^ F ^ a *K 2 (F * > a *) 9 "^ F *> a * ) ) 
V loo) V > ) £^™ K> iW 3\ >n 3 \ , A -r Z-j^ vv f tj (F*,a*) f tj (F*,a*) ^(F*,a*) ' 

(3.12) - g ^ ¥1^ ^/)) 

7,(5=1 ' ' 

with fij(F,a) = a 2 C t (F l )C J (F J )^^-, ^(F,a) = yftKV and A"" = [d^d 2 ]-\ 

Note that as opposed to other asymptotic methods presented in the literature, this formula is exact 
at t — 0. The zero-order formula (independent of the time t for Gi{t), Pij(t) constant) was derived 
for a general multi-dimensional local volatility model by 01 Moreover, in the expansion, we 
assumed that time is small but we made no assumption that Fk is close to the spot libor or that 
the volatility of volatility is small. 

3.4. Asymptotic Smile. The asymptotic smile can be derived in two steps from the asymptotic 
local volatility: first, we have 



(3.13) ds= y> ^(»o,w 



- Km) 

;oc( s o,i) 



and doing a change of local time t' = J Q cr"^.(so, u) 2 du, we now obtain the associated local volatility 
model for the swap rate 

(3.14) ds aP = al c {s,t)dWl 



(i 
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Secondly, we know that there is a one-to-one correspondence between this local volatility and the 
smile ^21 given at the first-order by (C(/) = ^if c (s, t = 0)) 



a/3 
T BS 



(K,T a 




(3.15) +12 



T 

C 3 (/a„) 



J«o C(/') 



"(I 



24 



' — w ' — ( — ^ at ^ \ 

C(/au) C{fav) 



0) 



)l/ a „ = ifl+£) 



4. SABR-LMM model 



We have seen that the asymptotic local and implied volatilities can be computed if we know the 
geodesic distance and a paramctrization of geodesic curves on M. n+1 . This is the case for the 
hyperbolic space H n for all n. This manifold has a lot of important properties. As such, it appears 
to be the perfect toy model (usually its Lorcntzian version AdS/dS) in a number of domain: 
chaos, cosmology, string theory, .... In the first part, we present our BGM-LLM-SABR model and 
show that the underlying geometry is H ,l+1 (with n the number of forward Libor rates). Using 
this connection, we will find an asymptotic local volatility and an asymptotic swaption implied 
volatility. 



4.1. Dynamics. We introduce the model SABR-LMM, given by the following SDE under the 
spot Libor measure Q (associated to the numeraire Bd(t) = nf=i *(! + T jFj(Tj — l))P(t,Tp^_i) 
where (3{t) = m if T m ^ 2 < t < T m _i) 



(4.1a) dF k = a 2 B k {F,t)dt + o k {t)aC k {F k )dZ k 

(4.1b) da = isadZ n+ i] , dZidZj = pij(t)dt i,j = 1, • ■ ■ ,n+ 1 

with 

(4.2) C k (F k ) = </> fe if fc 



(4.3) B k (F,t) = Yl 



The functions C k {F k ) have been scaled by <p k and therefore we can impose that a k (0) = 0. The 
stochastic equation for a was written in the spot Libor measure in order to get a SDE of a specific 
underlying swap s a p or a forward bond. Under the forward swap measure Q"' 3 , we have 



(4.4a) dF k = a 2 b k (F,t)dt + a k (t)aC k (F k )dZ k 

(4.4b) da = -v 2 a 2 b a (F, t) + vadZ n+1 \ , dZidZj = pij (t)dt i, j = 1, • • • , n + 1 
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BGM parameters 


a, b, c, d, p L , Sa, 5 b 


Cev parameters 


Pi, i = 1, ••• ,n 


SABR parameters 


a, v, p ia i = 1,- • • , n, 



Table 2. SABR-LMM: 9 + 3n parameters 



with 

(4.5) b k (F,t) 



J=«+l ^ v y i=min(t+lj+l) 



(l+T-i-Fi) 



2= Qf+1 

/3 



(4.6) C a0 (t) = Yl T iP(t,Ti) 

(4.7) b a (F,t) 

and with Wj(t) : 



+ r k F k (t) 



j2 T Mt)i2— Fh)Pka(Xk{t) 

=a + l k=l 



Uk=l (l + r,F k , 



2^ j = a + l 11 fc=l (l + r fc F fc ) 



Note that the forward-rate dynamics under the forward measure Q k is much simpler and given by 
the following stochastic differential equations (SDE) 



(4.8) 



dF k (t) = a k {t)aC k {F k )dW k , dW k dW p = p kp (t)dt 



As it is the case for the BGM model, we can use a piecewise parametric form or a functional form 
for the serial volatilities ui (t) and the correlation pij (t) (here full rank) as 



(4.9) 
(4.10) 



ai(t) = NiKaiTi-! -t) + d)e- b(T - 1 - t) + c] Vf < T;_i 

Pij(t = 0) = p L + (1 - pL)e~ {&A - SB minlT^T^m^-T^ 



The constants Ni are fixed such as eri(0) = 0. The model depends on 9 + 3n parameters (see Tab. 
2) which are calibrated on the swaption matrix. In the next subsection, we derive the metric, the 
geodesic distance and the Abelian connection underlying this model. 



4.2. Hyperbolic geometry. By definition, the infinitesimal distance (at t = 0) between the point 
x a and x a + dx a (l3~4alt (ds 2 = ££/i)=i g a pdx a dx ) is given by (p lj = [p'% , = (1, • • • ,n) 
and p m = are the components of the inverse of the correlation matrix p) 

After some algebraic manipulations, we show that in the new coordinates [x k ] k= i... n+ i (L is the 
Cholesky decomposition of the (reduced) correlation matrix: \p\i t j=x— n = [LL 

(4-12) x k = Y. vikl \ 7^7T+E^ a ' k = l,-- ,n 

4=1 Jf ? Ci ^i> i=l 

n 

(4.13) x n+1 = (l-^p«V%)5 fl 



8 



PIERRE HENRY-LABORDERE 



the metric becomes 

(414) ds 2 2 (! - Eh pV^O EIU ^ + M±i 

Written in the coordinates [xi], the metric is therefore the standard hyperbolic metric on H™ +1 

(modulo a constant factor 2 ^ ^'• 3 - k =£ p p ^- ), In order to compute our saddle-point we 
need the geodesic distance which is given by |19) 

Proposition 4.2.1. The geodesic distance d(x,x') on H™ +1 is given by 
(4.15) d(x,x°) = cosh- x [l + ^ l = x yXl ~ ° ] 

Using the geodesic distance on H™ +1 between the points x = ({F}k,a) and the initial point 
x° = {{F°} k , a) {q l = Jpo c7§7) ) g ivcn b y 

M o, ^(i-s^vv-vq)* ,_ lri , ^e"^i^^+ 2 ^; =1 p j %- + (^«) 2 i 

d(X,X ) = COsh H - ; == — ; r 

we derive the following non- linear equations l|3.11|l satisfied by the saddle-point a*(s), q*(s) which 
implicitly depends on s, the swaption strike: 



1 



»=1 ij'=l 



(446a) a*(s) 
(446b) 

(2^ a +2^^) d(aMg '* }) = -Aa(l - Y p pa P qa Ppq )^\* 



with 

(447) awfc(d(a*, {?*}) = 1 + 



(l-Ey=i/'V'%)a 
and with 

p*i-/3i _ pot-Pi 

(4-18) q* = <t>7\ 1 * 

(449) g* = C'MS) , A = 1 

The saddle-point is determined by solving these non-linear equations (|446() and an approximation 
(which could be used as a guess solution in a numerical optimization routine) is found by linearizing 
these equations around the spot Libor rates (i.e. qi = 0) 



-2v 2 (s- s ) - 2v V" Pij p> a 
(4.20a) A*(s)a 2 = — ^ ° qa ^ : + °((« - «o) a ) 

(1 - £p,g=l P Pa p qa p P q) 22iJ=l Pij^iUj 

(4.20b) = E |; lP ^ (S ' S0) + o((s - so) 2 ) 

with = (gj = 0). Note that when the strike is close to at-the- money, the saddle-points are 
close to the spot Libors and a* = a. 
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Moreover, by using the explicit expression for the hyperbolic distance, the Van-Vlcck-Morette 
determinant is 

(4.21) A(F, a, a) = - 

^/cosh 2 {d( ai F\a))-l 



4.3. Connection. The Abelian connection is given by l|3.4bj) 1 



(4-24, A, = ^g^^.afiffil).^,,,), 

(4.25) A. = i[^(gM_S2dffi ) _„ p .. r(i r, )] 
where we have used that 

(4.26) ^ = 2 ^ det ^ 



^niu^) 

Finally, the Abelian 1-form connection is 

A = \ D(7^t - ^P^)(" E ^ + - ^ E ^ + ^ QQda ) 

j=i 3V j'y i=1 ,, =1 

In order to compute the log of the parallel gauge transport ln('P)(a, q\a) = J c A, we need to know 
a parametrization of the geodesic curve on H ,l+1 . However, we can directly find ln(V)(a, q\a) if we 
approximate the drifts b k (F,t) by their values at the Libor spots (and t = 0). A similar approxi- 
mation was done in the Hagan-al formula as was shown in |12) . Modulo this approximation, 

i " hi(F° n\ d C (F°) " " 
ln(V)(a,q\a) ~ ^[(-J'i _ 3 ^ J )(^ P tf « + ~ "))] - 6 a (F°, 0)(^ + p QQ ( 

j=i i ' »=i i=i 

4.4. Asymptotic Smile-Summary. The asymptotic local volatility is given by H3.12[) 

d^fa^F) , 9 _^/ y (a,F)a^(a,F) 



«£) 2 M = E PiMt)°i(t)f ij (a,F)(l + 2T J2 A ^(( f ( F) 

i,j=l ^,f=l ' ' 

with (a* 2 (s), {f *}i(s)) the saddle-point satisfying the equations (|4.16|l and approximated by l|4.20[) 
and 



l 

(4.22) A* = ^{V -hCidiCi) 

u 2 a?b a (F,t) 

(4.23) A a = 
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fij(a,F) 
d(a,F) 
ln(T')(a, q\a) 
A(F, a, a) 



dR dFi 



cosh 1 \1 



v 2 E?,=i P»W + ^ ELi P 3a <lj + (<* - ") 2 



2(1 -££,-=! pV a /»«)< 



— )("L' J * 

j=1 JV J ; i=l 

d(o,F|a) 



; = 1 



v'cos/i^a, F|a)) - 1 
2T i cfe£[ / o]-2 



Note that this expression is exact when r goes to zero. The smile at the first-order is then obtained 
by plugging the above expression into i|3.15|l . 

Remark 4.4.1 (Libor CEV model). Note that our model for v goes to zero (and a = 1) gives the 
Andersen-Andreasen CEV libor model (with different CEV parameters for each libors) and the 
above expressions degenerates into 

, ds a/3 ds a/3 



(4.27) 
(4.28) 

(4.29) 

(4.30) 
(4.31) 



fij(F) 
d(F) 

ln{V){q) 
A(F,F°) 



en-; >([,{!■;,) 



dF t dF, 



2 P ij 1iti 

^ V(F°,0) d 3 C 3 {F») " 

^ [ Ci(F?) 2 J ^ 

j=l JV 3 ' i=l 

1 

2^det[pY^- 



P lj qi 



with the saddle-points (|4.1(i|) satisfying the non-linear equations (modulo the constraint s a p = s) 
(4.32) 



4.5. Comments and Numerical Tests. It is interesting to note that for n = 1, i.e. for a caplet, 
the caplet asymptotic smile reduces to the classical SABR formula by construction. 

Moreover, the asymptotic local volatility is given at the zero-order by 

Kf) 2 (M)= E ^(^(t)a,(i)a* 2 (F*)a(F;)a(F;)^(^)^(F*) 



with F* depending implicitly on s via l|4.16[) . At this stage, it is useful to recall how a similar 
asymptotic local volatility is derived using the "freezing" argument. The forward swap rate satisfies 
the following SDE in the forward swap numeraire Q Q/3 



(4.33) 



k=l 



ds a p = y^ j ^^-Uk{t)aC k (F k )dZ k 
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The "freezing" argument consists in assuming that the terms d gp® an d c(F ) are a l mos t constant. 
Therefore, the SDE (|4.33|) can be approximated by 



(4-34) ds aP = g ^ {F °)a k (t)a90±C k (s)dZ k 

and the local volatility is 

We can reproduce this formula for the swaption smile at-the-money 2 as the saddle-point Libor 
rates coincides with the spot rates. This is not the case for in/out-the-money swaption. Therefore 
our expression (exact at the zero-order) shows that the freezing argument is no longer correct 
when we try to fit a swaption implied smile in/out-the-money. In the following, we have tested 
our asymptotic swaption formula at the zero-order (Formula Fl) against the Andersen- Andreasen 
asymptotic formula (Formula F2) [Q in the case v = 0. The accuracy of these approximations 
are examined using Monte-Carlo (MC) prices as a benchmark. Following |16|. we consider five 
scenarii (see Tables 3-4-5-6-7). In the following tables, the implied volatility is reported and the 
numbers in brackets are the errors (in basis points i.e. true volatility times 10 4 ) corresponding to 
the implied volatility computed using the Fl or F2 formula minus the MC implied volatility. An 
x x y swaption has an option maturity of x years, a swap length of y years and a tenor of one year. 
We set a time-step for Monte-Carlo S = 0.125 and 2 16 paths 3 . Our formula Fl is more accurate 
than F2. 



Swaption 


strike 


MC 


Fl 


F2 




4%(ITM) 


22.42% 


22.41% (-1) 


22.61% (19) 


5 x 15 


6% (ATM) 


20.33% 


20.41% (8) 


20.46% (13) 




8%(OTM) 


18.92% 


18.93% (1) 


19.01% (10) 




4%(ITM) 


22.41% 


22.51% (11) 


22.67% (26) 


10 x 10 


6% (ATM) 


20.38% 


20.41% (3) 


20.50% (12) 




8%(OTM) 


18.93% 


18.93% (-1) 


19.05% (12) 



Table 3. Scenario A: Libor volatility A 4 (t) = 5%. Libor ^(0) = 6%. = 0.5 



Swaption 


Strike 


MC 


Fl 


F2 




5.08%(ITM) 


18.12% 


18.20% (8) 


18.17% (5) 


5 x 15 


7.26%(ATM) 


16.51% 


16.61% (10) 


16.63% (12) 




9.44%(OTM) 


15.38% 


15.38% (0) 


15.56% (18) 




5.55%(ITM) 


17.80% 


17.81% (1) 


17.89% (9) 


10 x 10 


7.93%(ATM) 


16.26% 


16.33% (7) 


16.38% (11) 




10.31%(OTM) 


15.17% 


15.19% (2) 


15.32% (15) 



Table 4. Scenario B: Libor volatility &U = 0.25(0.17 + 0.002(T^i - t))L?dW. 
Libor Li(0) = log{a + bi), L (0) = 5%, Li 9 (0) = 9% . = 0.5 



An at-the-money swaption (ATM) has a strike K equal to the spot rate s Q( 3(0) and an out-of-the money (OTM) 
(resp. in-the-money (ITM)) swaption has K < s a p(0) (resp. K > s a p(0)). 
3 We have used a predictor-corrector scheme with a Brownian bridge. 
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Swaption 


Strike 


MC 


Fl 


F2 




5.08%(ITM) 


14.89% 


14.97% (8) 


15.08% (19) 


5 x 15 


7.26%(ATM) 


13.73% 


13.79% (4) 


13.81% (8) 




9.44%(OTM) 


12.92% 


12.91% (-1) 


12.92% (0) 




5.55%(ITM) 


14.52% 


14.53% (1) 


14.64% (12) 


10 x 10 


7.93%(ATM) 


13.33% 


13.38% (5) 


13.40% (7) 




10.31%(OTM) 


12.51% 


12.51% (0) 


12.54% (3) 



Table 5. Scenario C: Libor volatility dL l = 0.25(0.17 - 0.002(Ti_i - t))L?dW. 
Libor Li(0) = log{a + bi). L (0) = 5%, L 19 {0) = 9%. = 0.5 



Swaption 


Strike 


MC 


Fl 


F2 




5.08%(ITM) 


19.19% 


19.33% (14) 


19.38% (19) 


5 x 15 


7.26%(ATM) 


17.59% 


17.72% (13) 


17.75% (16) 




9.44%(OTM) 


16.46% 


16.49% (3) 


16.61% (15) 




5.55%(ITM) 


18.92% 


18.94% (2) 


19.06% (14) 


10 x 10 


7.93%(ATM) 


17.31% 


17.39% (8) 


17.45% (14) 




10.31%(OTM) 


16.18% 


16.21%(3) 


16.32% (14) 



Table 6. Scenario D: dL; 



0.05Lf( . b ' l(t) dW l 

1 V b «(*) 2 +M*) 2 



-dWo 



— ; = n, v v 2 I • 

Vfc-l(t) 2 +fc.2(t) 2 *> 

b a (t) = ^e-*i( T «-i-*) + e e -*am-i-t) j bi2 (t) = V /T— ^(Jl-i-O, p = 0.99, 
6 = -0.99, ki = k 2 = 0.54. Libor L ; (0) = log(a + bi). L o (0) = 5%, Li 9 (0) = 9%. 
(3 = 0.5 



Swaption 


Strike 


MC 


Fl 


F2 




5.08%(ITM) 


33.09% 


33.65% (56) 


34.24%(115) 


5 x 15 


7.26%(ATM) 


29.47% 


29.96% (49) 


30.23%(76) 




9.44%(OTM) 


26.92% 


27.14% (22) 


27.49% (57) 




5.55%(ITM) 


31.75% 


32.41% (66) 


33.33% (158) 


10 x 10 


7.93%(ITM) 


28.47% 


28.88% (41) 


29.37%(91) 




10.31%(OTM) 


26.01% 


26.18% (17) 


26.68% (67) 



Table 7. Scenario E: Scenario D with (3 = 0.3 



In ^2], w e explained how to derive a general asymptotic smile for any stochastic volatility model 
using this geometric framework. As an application, we derived an asymptotic smile for a SABR 
model with a mean-reversion term. In the following, we try to consider some natural extensions 
of our SABR-BGM model where we add a non trivial drift to the volatility process. The only 
modification comes from the expression of the Abelian connection. 



4.6. Extensions. Under the spot Libor measure, we assumed that the volatility follows the process 
(4.35) da = -va 2 4> a (a)dt + vadZ n+1 
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with ip a (a) a general analytical function of a (the scaling a 2 in front of ij) a (a) has been put for 
convenience). After some algebraic computations, we derive the new Abelian 1-form connection 

A = I Ek^SS - djCi 2 Fi) )("itp ijd * + p aJda ^ - ( 6 °( F ' *) + ^(°))("E piadqi + paada) 

j=i i\ i> i=l i=l 

Using a similar approximation as before, i.e. Cj(Fj) ~ Cj(Fj) and ip a (a) ~ -0 a (a), we obtain for 
the parallel gauge transport 

i " K>(F° 0\ d-C (F°) " 
Zn^MM ~ £[(_J'J - J ^ j; )(^^^+p^(a - a))] - (6 a (F°,0) 

j=l ^ J' i=l 
n 

(4.36) + ^(a))(i/£y° ft +//»°(a-a)) 

i=l 

Finally, the smile is obtained using our general formula 1)3. 15|) . Note that the metric and the 
geodesic equations remain unchanged when we only modify the drift terms. 



5. Conclusion 



In this short note, we have introduced a LMM model coupled to a SABR stochastic volatility 
process. By using the heat kernel expansion technique in the short time limit, we have obtained 
an asymptotic swaption implied volatility at the first-order, compatible with the Hagan-al classical 
formula for caplets. Moreover, we have seen that this exact expression (when the expiry is very 
short) is incompatible with the analog expression obtained using the freezing argument. 
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